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1. Introduction 

The strong coupling constant and the quark masses constitute fundamental parameters of the 
Standard Model. It is an important task of lattice QCD to determine these parameters using inputs 
at low energy scales. In the course of evaluating these fundamental parameters we need the process 
of renormalization in some scheme. The MS scheme is one of the most popular schemes, and 
hence one would like to evaluate the running coupling constant and quark masses through input of 
low energy quantities on the lattice and convert them to the MS scheme. A practical difficulty in 
this process is called the window problem that the conversion should be performed at much higher 
energy than the QCD scale. At the same time the renormalization scale pL should be kept much 
smaller than the lattice cut-off to reduce lattice artifacts: l/L <C Aqcd <C /X <C I /a. 

The Schrodinger functional (SF) scheme ^ ||] is designed to solve the window problem. A 
unique renormalization scale is introduced through the box size L in the chiral limit and the scheme 
is mass independent. A wide range of renormalization scales can be covered by the step scaling 
function (SSF) technique. This matches our goal to obtain the coupling constant and quark masses 
in the MS scheme. The SF scheme has been applied for evaluation of the Nf = 2 + 1 QCD coupling 

For the quark mass renormalization factor the SF scheme has been applied for Nf = 2 QCD 
[Q]. At low energy scales of /J. ~ 500 MeV, where physical input is given, we expect the strange 
quark contribution to be important in addition to those of the up and down quarks. Thus the aim of 
this report is to go one step further and evaluate the quark mass renormalization factor in Nf ■ = 2 + 1 
QCD 1 . Our goal is to renormalize the bare light quark masses evaluated in two recent large-scale 
Nf = 2 + 1 lattice QCD simulations ^ [7[ [|[] and derive the renormalization group invariant (RGI) 
quark mass M. Once the RGI mass M, which is scheme independent, is evaluated, the conversion 
into the MS scheme can be carried out perturbatively. 

We therefore first derive the renormalization factor Zm(§o), which converts the bare PC AC 
mass at a bare coupling go to the RGI mass. Derivation of the renormalization factor proceeds in 
three steps [Q, g] 

Z A (g Q ,a/L) m{\/L n ) M 

^M{go,a L max — — . (1.1) 

Zp(g ,a/L max ) m(l/L max ) m(l/L n ) 

The first factor Z A /Zp renormalizes the bare PC AC mass in the SF scheme at an appropriately low 
energy scale L max . The second factor m(l/L,,)/m(l/L max ) represents the running of the renormal- 
ized mass from L max to a high energy scale L n , which is evaluated non-perturbatively in the SF 
scheme. The last factor M/m(l /L n ) represents the running from L n to infinitely high energy scale 
and is evaluated perturbatively for an appropriately high energy scale L n . 

2. Step scaling function 

We adopt the renormalization group improved Iwasaki gauge action and the non-perturbatively 
0(a) improved Wilson fermion action with the clover term. The Dirichlet boundary condition for 
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the spatial gauge link is set to 

U k (x)\ X0=0 = U k (x)\ X0=T = l (2.1) 
and the twisted periodic boundary condition of the fermion fields in the three spatial directions is 
y(x + Lk) = e w yr(x) , W(* + Lk) = e~ ie y(x) , 0=0.5. (2.2) 

We prepare seven renormalized coupling values to cover weak to strong coupling regions 
[g]. For each coupling we use three box sizes L/a = 4,6,8 to take the continuum limit. At the 
three lattice sizes the values of j3 and K were tuned to reproduce the same renormalized coupling 
keeping the PCAC mass to zero. On the same parameters we evaluate the pseudo scalar density 
renormalization factor Zp(go,a/L) and Zp(go,a/(2L)) at two box sizes. 

The renormalization factor Zp is defined in terms of two-point functions [Q] of pseudo scalar 
density at bulk and at boundary. Taking the ratio of two renormalization factors at two renormal- 
ization scales /i = 1/L and 1 / (2L) we get the step scaling function (SSF) on the lattice 



a\ _ Z P (g 0l a/(2L)) 

Li . 



' LJ Zp(g ,a/L) 



(2.3) 

g 2 (L)=«,m=0 



where mass independent scheme in the massless limit is adopted. 

We perform a perturbative improvement of the SSF before taking the continuum limit. For this 
purpose we need a perturbative evaluation of the lattice artifact in the SSF 

Op(u,a/L) = p-r . (2.4) 

Op{u) 

Instead of calculating 8p at one and two-loop level perturbatively we calculate SSF's directly by 
Monte-Carlo simulations at very weak coupling j8 > 10. The artifact is fitted to a polynomial form 
for each a/L, 

l + 8 P (u,a/L) = 1 + d y (a/L)u + d 2 {a/L)u 2 . (2.5) 

Since the quadratic fit provides a reasonable description of data [^] we opt to cancel the 0(a) 
contribution dividing out the SSF by the quadratic fit. 

Scaling behavior of the improved SSF is plotted in the left panel of Fig. [T[ Almost no scaling 
violation is found. We performed three types of continuum extrapolation: a constant extrapolation 
with the finest two (filled symbols) or all three data points (open symbols), or a linear extrapolation 
with all three data points (open circles), which are consistent with each other. The constant fit was 
employed with the finest two data points to find our continuum value. 

The RG running of the continuum SSF is plotted in the right panel of Fig. [jj A polynomial fit 
of the continuum SSF to third order yields 

Gp(u) = \+pou + piu 2 + p 2 u 3 , (2.6) 
p x = -0.002826, p 2 = 0.000031 (2.7) 

fixing the first coefficients po to its perturbative value. The fitting function is also plotted (solid 
line) together with the two loops perturbative running (dashed line). 
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Figure 1: The SSF on the lattice with its continuum extrapolation at each renormalization scale (left) and 
the RG flow of the SSF (right). 
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Figure 2: The non-perturbative running mass in the SF scheme. Solid line is a perturbative running with 
two and three loops RG function for the mass and the coupling. 



Multiplying the SSF according to a sequence of couplings w,-, which differ by a factor two in 
the renormalization scale, starting from uq = g 2 {L max ) we get the non-perturbative running of mass 

" . . m(l/L n ) 

In Fig. H| we plot the non-perturbative running mass m(\/L n )/M in units of the RGI mass as a 
function of the scale /x/Asf, where jj. = l/L n . 



3. Zp and Za at low energy scale 

We evaluate the renormalization factor Zp at the same bare coupling /3, 1.83, 1.90 and 2.05 
adopted in the large scale simulations [^, ^ §J. The reference scale L max is given by the box size 
we used in this evaluation. The renormalized coupling g 2 (L max ) should not exceed our maximal 
value 5.13 for the coupling SSF significantly. We use the box size of L/a = 4 for j8 = 1.83 and 
1.90 to define L max and L/a = 6 for j8 = 2.05. We adopt the lattice spacing a as an intermediate 
scale. 
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Figure 3: j3 dependence of Za(#o) with and without disconnected diagrams (left) and comparison with 
perturbative results (right). 



The remaining ingredient of the renormalization factor ( LI ) at low energy is Za of the axial 
vector current. We calculate the renormalization factor according to the procedure in Ref. [^] 
through the axial Ward-Takahashi identity, which is applicable to small non-vanishing PCAC mass. 
The box size L for the axial current need not coincide with L max since Za is scale independent in 
the continuum limit. 

In this paper we adopted the following box size at each /3 to define Za'. 12 3 x 30, 8 = 0.5 
at J8 = 1.83, 10 3 x 24, = at J8 = 1.90 and 12 3 x 30, 6 = 0.5 for j3 = 2.05. The disconnected 
contribution ^ is included. In the left panel of Fig. || we plot the values for Za according to our 
definition (filled circles), together with those without the disconnected contribution (open symbols). 
Scattering of points starting around j8 = 1.95 indicates that lattice artifacts are increasingly large 
in our data for large lattice spacings. In the right panel of Fig. || we compare our Za(§o) with 
perturbative behavior (solid line) and results from the tadpole improved perturbation theory (star 
symbols). 



4. RGI mass renormalization factor 



We derive the renormalization factor Zm for the RGI mass, which is intended to renormalize 
the bare PCAC masses obtained in the two large scale simulations. The three hadron masses, 
m K , nix, tfiQ were used in Ref. [[7|, |8|] to determine the light quark masses and the lattice spacing. 
Two choices m n , m p mk or m K , m p m§ were adopted in Ref. [Qj. The results for Zm are listed in 
Table [|. Also listed are the renormalization factors Z^ s (jS , jU =2 GeV) in the MS scheme at a 
renormalization scale \i = 2 GeV. We emphasize that the renormalization factor here is defined in 
terms of the renormalization group functions for three flavors. 

The error in the renormalization factor includes all the statistical and systematic ones except 
for that from the choice of the reference scale L max . We tried a rough estimate of 0((a/L max ) 2 ) 
error and we find a few percent effect at /3 = 2.05, while the magnitude may increase to a 10 % 
level at lower values of j8 . However, a firmer conclusion requires the step scaling function at the 
couplings stronger than those explored in the present work. 

As the last step we apply our renormalization factor to the bare PCAC masses in Refs. [||, [7J, ||] 
to obtain values for the renormalized quark masses. We notice that the chiral perturbation theory 
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p 


Zm 




PT(tad) Z™ s (K) 


(0) 


PT(tad) Z™ s (0) 


1.83 [g] 


1.33(15) 


1.04(12) 


1.07161 


1.04(12) 


1.07019 


1.90 [g 


1.693(46) 


1.315(35) 


1.09973 


1.315(35) 


1.09955 


2.05 [| 


1.862(41) 


1.417(29) 


1.14487 


1.416(29) 


1.1446 


1.90 g 


1.693(46) 


1.344(36) 


1.11322 






1.90 [g 


1.693(46) 


1.347(36) 


1.11322 







Table 1: Zm for the RGI mass and Z™ s {2 GeV) in the MS scheme. (K) or (<j>) means which meson mass is 
used for physical scale input. Perturbative results with tadpole improvement are also listed. 

was used to extract the physical quark mass in Ref. H 0, while the reweighting technique was 
adopted in Ref. ^ to evaluate the quark mass directly on the physical point. The numerical results 
are given in Table |2| and are plotted in Fig. |] both for the averaged up and down quark mass (left 
panel) and for the strange quark mass (right panel). For the old CP-PACS/JLQCD work of Ref. [^] 
results with K and <p meson input are plotted (filled squares and circles) together with perturbatively 
renormalized masses using tadpole improved renormalization factor (open squares and circles). 
The upward triangle represents the result for the more recent work of PACS-CS evaluated directly 
on the physical point using the reweighting technique [j|]. The downward triangle represent a result 
of chiral extrapolation from the simulation point reaching down to m n = 155 MeV [0]. 



p 






Mf GI 


< s 


1.83 [§] 
1.90 
2.05 [§] 


3.30(38) 
4.47(17) 
5.29(24) 


2.59(30) 
3.47(13) 
4.02(18) 


86(10) 
115.3(4.4) 
136.1(6.5) 


67.7(7.9) 
89.5(3.4) 
103.6(4.9) 


1.90 0] 
1.90 § 


3.84(16) 
3.49(34) 


3.05(12) 
2.78(27) 


110.5(4.0) 
109.0(3.0) 


87.7(3.1) 
86.7(2.3) 



Table 2: Non-perturbatively renormalized mass for the averaged up and down quark and for the strange 
quark. M RGI is the RGI mass and m MS is that in the MS scheme at a scale fi =2 GeV. The unit is in MeV. 
/f-meson mass is used for physical input. 

It is disappointing that the old CP-PACS/JLQCD results do not exhibit a better scaling behavior 
by going from perturbative to non-perturbative renormalization factor. However, we should note a 
significant change in the average up and down quark mass with the recent PACS-CS work (filled tri- 
angles). This represents a systematic error due to chiral extrapolation of the old CP-PACS/JLQCD 
work whose pion mass reached only m K ~ 500 MeV. We should also note that the renormalization 
factor Za involves a large uncertainty at j8 = 1.83 which is not reflected in the error bar of Fig. |]. 
We feel that results at j8 = 2.05 using simulation with physical pion mass are needed to find the 
convincing values for light quark masses with our approach. 

5. Conclusion 

We have presented a calculation of the quark mass renormalization factor for the Nf = 2 + 1 
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Figure 4: Scaling behavior of (left) and m^ (right). 



QCD in the mass independent Schrodinger functional scheme in the chiral limit. We calculate the 
SSF of the running mass on the lattice and applied a perturbative improvement. We find that the 
step scaling function shows a good scaling behavior and the continuum limit may be taken safely 
with a constant extrapolation. The non-perturbative SSF turned out to be almost consistent with 
the perturbative two loops result. We then derive the renormalization factor of the pseudo scalar 
density at low energy and that of axial vector current. Multiplying these factors we evaluate the 
non-perturbative renormalization factors for the RGI mass and that in the MS scheme at a scale 
jU = 2 GeV. 

Applying our non-perturbative renormalization factor to the present PACS-CS simulation re- 
sult evaluated directly at the physical point yields m^{pL = 2 GeV) = 2.78(27) MeV for the aver- 
age up and down quark, and m^ s (/i = 2 GeV) = 86.7(2.3) MeV for the strange quark. Simulations 
at weaker couplings under way will tell if these values stay toward the continuum limit. 

This work is supported in part by Grants-in-Aid of the Ministry of Education (Nos. 10143538, 
20105001, 20105002, 20105003, 20105005, 20340047, 20540248, 20740139, 21340049, 22105501, 
22244018, 22740138). 
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